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Abstract 

We investigate zero-temperature dynamics on the hexagonal lattice EI for the 
homogeneous ferromagnetic Ising model with zero external magnetic field and a 
disordered ferromagnetic Ising model with a positive external magnetic field h. We 
consider both continuous time (asynchronous) processes and, in the homogeneous 
case, also discrete time synchronous dynamics (i.e., a deterministic cellular au- 
tomaton), alternating between two sublattices of H. The state space consists of 
assignments of —1 or +1 to each site of H, and the processes are zero-temperature 
limits of stochastic Ising ferromagnets with Glauber dynamics and a random (i.i.d. 
Bernoulli) spin configuration at time 0. We study the speed of convergence of the 
configuration o l at time t to its limit <7°° and related issues. 
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1 Introduction 

In this paper, we consider a number of continuous time Markov processes cr*, t > 0, with 
state space S = { — 1, +1} L consisting of assignments of —1 or +1 to a regular lattice L. 
In general, we take L to be H, the hexagonal lattice in the plane, but occasionally we refer 
to results on Z d or other lattices. Later on, we will also consider a related discrete time 
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process on EL A state a £ 5 will also be called a spin configuration. All the continuous 
time processes we consider are nearest-neighbor interacting particle systems in which the 
dynamics may be constructed by means of independent, rate 1, Poisson "clock" processes 
at the sites x of L. If the Poisson clock at x rings at time t, then a spin flip (crj. + = — cr* — ) 
is considered; whether or not a flip actually occurs depends on the values of a l ~ and <r*~ 
for each neighbor y of x. 

For the simplest model, the homogeneous ferromagnet, there is a flip with probability 
1 (or 1/2 or 0), if at t— , o~ x agrees with less than half (or exactly half or more than 
half) of its neighbors. Note that the probability 1/2 case does not occur in lattices like 
EI where every site has an odd number of neighbors so that there cannot be a tie among 
the neighbors. We denote by u a realization of the clocks rings (and tie-breaking coin 
tosses, if needed) and by P u the corresponding probability measure. In all our models 
we choose cr° according to a probability measure P a a corresponding to i.i.d. cr°'s with 
P a o(a® = +1) = A G [0,1]. In the Ising model context, A = 1/2 is the most important 
special case and corresponds to an initial "quench from infinite temperature" . 

In the physics literature on zero-temperature dynamics for the homogeneous ferro- 
magnet with A = 1/2, an important quantity is the P a o x P u probability p(t) that the 
origin has not flipped at all by time t. In this situation, it has been known for a long time 
that on Z, almost surely (a.s.) every site flips infinitely often and hence p(t) — > as 
t — > +oo. More recently it was shown |], |7j that on Z, p(t) ~ t -3 / 8 and numerically || 
that on Z 2 , p(t) ~ t' 6 ^ with 0(2) « 0.22. This is consistent with the theorem [gTJ that 
also for d = 2, a.s. every site flips infinitely often. There is as yet no rigorous result for 
Z d with d > 2, but numerical evidence J3U] suggests that for d = 3, pit) ~ t~ e ^ with 
< 9(3) < +oo and raises at least the possibility that pit) — > p(oo) > for d > 4. If the 
latter occurs, it may still be that p(t) — p(oo) ~ t~ e ^ with < 9(d) < +oo. (We discuss 
the case A ^ 1/2 below.) 

The situation for lattices like H with no ties among neighbors is quite different. Here 
it is easy to see that p(oo) > (e.g., on EI, consider elementary hexagons with constant 
spin value at t — ) and a general theorem [2l| implies that a.s. every site flips only 
finitely many times so that cr°° = lim^oo a 1 exists. A natural question is then whether 
p(t) — p(oo) or its upper bound 



p(t) = (P a o x PJ) (origin flips after time t) 



tends to zero like t~^ or exponentially fast or something in between. 

Based on results for simpler quasi-one- dimensional lattices like Z x {0, 1} 
results for the homogeneous tree of degree three fll8| , it was speculated i] 
convergence to the fixed, i.e., absorbing, state <t°° may be exponentially fast: 



| and on 
31 that on H 



p(t) < exp{-(ct)}. 



(2) 



This possibility is also supported by numerical evidence of Howard [0, [Tj|. Although 
this was originally considered for A = 1/2, it seems reasonable that such a conjecture 
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should be valid for all A. One of the main results of this paper is a proof of the following 
stretched exponential estimate with exponent 7 = 1/2 for A not too close to 1/2 — see 
Theorems [1] and in Section 2 below: 

p(t) < exp{-(cT)}. (3) 

In this regard it is interesting to note that on Z d with d > 2, and A sufficiently close to 1 
(respectively, 0), it has been proved (12| that a 1 — > <j°° = +1 (resp., = —1) with a similar 



stretched exponential upper bound, but with 7 = j(d) 6 (0, 1), and a similar lower bound 
for d = 2. 

Closely related to bounds like (|2|) or @ is the spatially localized mechanism for fixa- 
tion, namely a x can no longer flip once it belongs to either a loop or "barbell" of constant 
sign in H. By loop we will always mean a simple loop (with no subloops). A barbell 
consists of two disjoint loops connected by a path and we regard a loop as a degenerate 
barbell. By studying the percolation properties of the final state a°° on the infinite lattice 
H, it follows - see Prop. |2.2| - that for any A, every site is in some <7°°-barbell of constant 
sign. Let R denote the diameter of the smallest such barbell. We conjecture, that for any 
A, 



(P ao x P^R > r) < exp(- C "r) (4) 

for some c" G (0, 00). For A not close to 1/2, this follows from the proofs of Theorems 
[1] and below. There would be an interesting application of the combination of (|3]) and 
(|j) concerning overlap distributions (cf. Sec. IX of p2|), which we will briefly describe at 



the end of this section after we introduce zero-temperature dynamics for disordered Ising 
models. 

Our results for the continuous time dynamics for the homogeneous ferromagnet on H 
are in Section 2. Then in Section 3, we analyze a discrete time dynamics on EI (alternating 
between and synchronous within each of two sublattices) that is the zero-temperature 
case of dynamics studied by Domany B and that has also been studied numerically by 



Nienhuis ||26|| . This is really a deterministic cellular automaton and we obtain much 
stronger results for it, such as exponential decay of p for any A, than we have obtained for 
the usual continuous time dynamics. We remark that some interesting results are obtained 
in m on the continuum scaling limit of the dependent percolation models generated by 
this cellular automaton when A = 1/2. Finally, in Section 4, we study continuous time 
dynamics for disordered Ising models on HI. 

In disordered Ising models, the flip probability is determined by a (formal) Hamiltonian 

H = - } u J x , y a x a y - h^a x , (5) 

(x,y) x 

where Yl{xy) denotes the sum over all pairs of nearest neighbor sites. We write Af(x) for 
the set of nearest neighbors of x. The homogeneous ferromagnet corresponds to the case 
where 
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• J x y = 1 for all pairs of nearest neighbors x, y, 

and we also take h = in that case, while disordered models correspond to the case where 

• the J Xt y's are i.i.d. random variables. 

We denote by Pj the distribution of coupling realizations J . When the Poisson clock at 
x rings, if 

AH x (a) = 2 22 Jx,yVx(Ty + 2hcr x (6) 

y&jV(x) 

is negative, then the flip is done with probability 1; if it is positive, then the flip is 
done with probability 0. If the energy change were zero, the flip would be done with 
probability 1/2; however for the disordered models, we will generally consider couplings 
whose common distribution is continuous, so that there is zero probability for zero energy 
flips to occur. For such disordered models with continuous distributions of couplings, it is 
also a consequence of a general result of that cr°° = lim^oo a 1 exists. Again, we will 
be interested in the t — > oo behavior of p(t), which here is the P a o x P w x Pj probability 
that the origin flips after time t. 

Similar stochastic processes on different types of lattices have been studied in various 
papers. See, for example, ||, 12, |14], ^T], |23|, 25| for models on Z d and |L8| for a 



homogeneous ferromagnetic model on T3, the homogeneous tree of degree three. Such 
models are also discussed extensively in the physics literature, including homogeneous 
and disordered ferromagnets as well as spin glasses, usually on the 7L d lattice. 

Exponential decay of pit) has been obtained by Newman and Stein |23) for the homo- 
geneous ferromagnet on the "ladder" lattice (Z x {0, 1}) and for continuously disordered 
models on Z. In both settings, the exponential decay is a result of a "blocking" condition 
present at time zero and due to the initial random spin configuration in the case of the 
ladder and to the structure of the random couplings in the case of the disordered ferro- 
magnets on Z. This blocking conditions break up the system into connected components 
that can have no influence on each other, making the dynamics fundamentally local in 
nature. The situation of the homogeneous ferromagnet on T3 is substantially different. 



It is shown in [[Tj| that, if the density of +1 spins at time is large enough, the system's 
agreement-inducing dynamics produces enough "fixated spins" to break up the system 
and ensure exponential decay of p(t). The blocking mechanism in this case is dynamical. 

We now turn to a discussion of overlap distributions for disordered (and homogeneous) 
models. Let us denote by fi°° the probability distribution of the random final state a°° 
(induced by P UQ x P u or by P ao xP w x Pj in disordered models) and let a 00 ' be a "replica" 
of ex 00 — i.e., let the pair (cr 00 ,^ 00 ') be distributed by x p,°°. Then, by the spatial 
ergodic theorem, it follows [E2| that the overlap random variable 



Q = lim - 

L-00 1 Ax, n 



E 



(7) 
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where denotes the L x L square centered at the origin and | • | denotes cardinality, 
is a constant g = g(X) (with g(l/2) = 0). It is however a priori possible (see [^] for a 



discussion) that for spin glass models, where the common distribution of the couplings 
is symmetric about zero, this triviality of the dynamical overlap distribution would not 
be so if one instead considered the limit (in distribution) as L — > oo of the finite volume 
overlap 



^ s |A7nir £ °?M°?M> (8) 

where o"°°[Ai] is the final state of the zero-temperature Markov process restricted to the 



finite volume Al (with some boundary condition). But it was conjectured in |22[ that 
Ql should rather converge in probability to the same constant g(X). To prove this would 
require that cr°° [A J — > a°° as L — > oo in such a way that a 00 [AJ and a°° only disagree in 
a small "boundary layer" (of area o(L 2 )). In Corollary |2.1| we do show that cr°°[A L ] — > 
and the extra uniformity needed would follow easily by combining conjectures and @. 
We have in fact this result for A not too close to 1/2, but the physically more interesting 
case is A = 1/2. Of course, the most interesting result would be for a spin glass rather 
than for a ferromagnet, homogeneous or otherwise. 



2 Homogeneous ferromagnet 

In this section, we consider the homogeneous ferromagnet on H with zero external mag- 
netic field, that is, J XjV = 1 for all pairs of nearest neighbors x, y and h = 0. Let A be the 
density of +1 spins in <j°. 

Our first result concerns the almost sure absence of percolation of both +1 and —1 
spins in the final configuration a°° (or in <j t ) when A = 1/2. We note that it has been 
proved |B| that in o~°°, the mean cluster size is infinite (see also Theorem^ below). 

Proposition 2.1. In the homogeneous ferromagnet, if A = 1/2 and t 6 [0, oo], then for 
almost every a and uj there is no percolation in o l of either +1 or —1 spins. 

Proof. First note that the measure // describing the state a 1 of the system at time 
t G [0, oo] is invariant and ergodic under any H-automorphism. This is so because the 
same is true for both P a o and and hence also for P a o x P w . Applying a result of Harris 



jl6| , p0| , we also have that // satisfies the FKG property, i.e., increasing functions of the 
spin variables are positively correlated (this follows from the FKG property of P a o and the 
attractiveness of the Markov process). Then it follows from a result of Gandolfi, Keane 



and Russo (stated in [O for Z 2 , but also valid for EF) that if percolation of, say, +1 sites 



were to occur, all the —1 clusters would have to be finite. Because of the symmetry of the 
model under a global spin flip, however, percolation of +1 sites with positive probability 
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implies the same for —1 sites. Then, using the ergodicity of the measure, we would see 
simultaneous percolation of both signs, thus obtaining a contradiction. □ 



Remark 2.1. Numerical evidence for the homogeneous ferromagnet ]7^/ suggests that 
there is plus (resp., minus) percolation in cr°° for all A > 1/2 (resp., A < 1/2). For such 
a result in the case of synchronous dynamics, see Prop. \3.1\ below. 



Remark 2.2. The proof of Proposition \2.i{ works in a more general context and was used 
in /fj/ to get a similar result for the homogeneous ferromagnet on 7? (see Proposition 3.2 of 
/d/j. The proof also shows that the symmetric Bernoulli product measure may be replaced 
by any distribution for cr° which (a) is symmetric under cr° — > — a , (b) is invariant and 
ergodic under any (nontrivial) H- automorphism, and (c) satisfies the FKG property. 

Let us call a configuration a G S locally fixated if for each x G H, there exists a (finite) 
lattice animal A x containing x such that for any b G S that coincides with o on A x , 
AH,y(a) > for every y G A x , i.e., y agrees with a strict majority of its neighbors. 

Proposition 2.2. In the homogeneous ferromagnet, for any X, for almost every a and 
uj, the final configuration cr°° is locally fixated. 

Proof. Let's first consider the case A = 1/2. In this case the claim follows immedi- 



ately from Proposition ^]T] for t = oo and the fact that a°° is a.s. fixated. In fact, for each 
x G H, it is enough to take A x to be the almost surely finite (e.g. plus) cluster at x. We 
now write cr°°(A) to indicate dependence on the parameter A and consider A > 1/2. Let's 
couple cr*(A) with cx*(l/2), where <r'(A) and a'(l/2) have the same dynamics realization uj 
and <7°(l/2) < cr°(A) where < indicates the natural partial order. From the attractiveness 
of the dynamics it follows that <7*(l/2) < cr*(A) for all t G [0, oo]. Hence, each site x such 
that (both cr£°(A) = +1 and) a x <3 (l/2) = +1 is contained in an almost surely finite cluster 
C' x with C' x C C x , where C x is the cluster at x in cr°°(A) and C' x is the cluster at x in 
o-°°(l/2). Then, we let A x = C' x . For a site x such that <r~(l/2) = -1, the cluster C' x 
that contains x in a™ (1/2) is almost surely finite and therefore surrounded by C x , the 
finite union of all neighboring plus clusters. (This union is actually a single (connected) 
cluster, as can be seen using arguments connected to those used in the proof of Theorem 
f§ but we will not use that fact.) Since a°°(l/2) < <x°°(A), C+ is still plus in a°°(A). We 
let A x = C' x U C x . This completes the proof for A > 1/2. The case A < 1/2 follows by 
symmetry. □ 

The property of a°° of being locally fixated means that the dynamics produces finite 
clusters that are stable for local reasons. This, in turn, has the following interesting 
consequence (see, for example, p2] for some physics motivation). Consider the sequence 
A L of squares and denote by cr*[Ai] the Markov process defined on the sub lattice A L nW 
(with some boundary condition). Then it is easy to see that for every x G H, and using 
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a natural coupling between the processes for different L's (see the proof of Corollary [2.1| 
below), the following is true almost surely: 

lim ( lim o*[A L ]) = lim o* = of. (9) 

This follows from a finite speed of propagation of information argument (an argument 



of this kind goes back to Harris [|Tj]] and is necessary even to prove that the stochastic 
process itself is well defined; for a reference, see |L7| or p. 119 of jiTi). It is not clear 
though whether the two limits (t — > oo and L — > oo) commute. In fact, on a homogeneous 
tree, interchanging the space and time limits can produce a different result (for example, 
choosing the plus boundary condition). But in our case, Proposition |2]2| implies the 
following corollary: 

Corollary 2.1. Almost surely, 

lim (lim oi[A L )) = lim (lim <t£[A l ]) = of. (10) 

L— >oo \t— >oo / t — >oc \L— >oo / 

Proof. For almost every initial configuration cr° and realization of the dynamics uj, 
cr°° (cr^, cu) is locally fixated and therefore, at time t = oo, x belongs to a lattice animal 
A x which is stable for local reasons. Call Ta x the time it takes for all the sites y G A x to 
fixate. < oo almost surely. Let us now introduce a coupling between cr'[A^] for every 
L and a 1 constructed as follows: 1) the dynamics realizations for cr*[A^] and a 1 coincide 
in A^nH, and 2) <r° = ct°[Al] Vx G Pi HI. If, for any given e > 0, one can take L large 
enough so that A x C A^ fl EI with probability at least 1 — e and so that moreover 

P^l A ;[^L]^a T / : \a G [A L ]=al L )<e, (11) 

where Ob is the configuration a G S restricted to the set B, the corollary follows. But for 
any e > 0, such an L = L(e) exists because of the finiteness of Ta x and the finite speed 



of propagation of information (once again, see |T7j or p. 119 of fL0|), and this concludes 
the proof. □ 



Before we state our next result, we need some notation. Call any minimal hexagon 
connecting six sites of HI a cell. Each site is contained in three cells. Given any connected 
subset A of HI, the internal energy of A (denote by e^) is defined as 

ZA = - Jx^x^y, (12) 

(x-y) 
x,y£A 

the sum being performed over pairs of neighboring sites, each pair counted once. 

Proposition 2.3. In the homogeneous ferromagnet, starting from any initial state (i.e., 
spin configuration) , each site flips at most 8 times. 
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Proof. By the translation invariance of the model, it is enough to prove the claim 
for the origin 0. Let A be the union of all sites belonging to the three cells that contain 
the origin. Then —15 < ca < 15. Every spin flip of the origin or of one of its three 
neighbors lowers by at least two units the internal energy of A, while a spin flip of any 
one of the other sites in A can lower the internal energy ca or leave it unchanged, but 
can never raise it. Therefore ca can never increase and moreover it strictly decreases by 
at least two units whenever the origin or one of its neighbors flips. Notice that after the 
origin has flipped the first time, every other spin flip of the origin has to be preceded by 
the spin flip of at least one of its neighbors. Thus every spin flip of the origin after the 
first one corresponds to a decrease in eA of at least four units. Given that the maximum 
change in ca (the energy when all edges are unsatisfied minus the energy when all edges 
are satisfied) is 30, the origin can flip at most 8 times. 

Remark 2.3. An analogous result holds true for the homogeneous ferromagnet on the 
ladder lattice (Z x {0,1}). The proof is the same once the corresponding cell has been 
properly defined. 

Recall that p(t) denotes the probability that the origin flips after time t. The following 
result says that, starting from a large enough density A of +1 sites, the time it takes for 
a 1 to converge to <t°° has a tail decreasing at least as fast as a stretched-exponential. We 
denote by p s c lte the critical value for independent site percolation on HI. The symbols c, c' 
and c" will denote generic positive constants, whose values may be different in different 
parts of the paper. 

Theorem 1. In the homogeneous ferromagnet, if X > pf te , then 

pit) < e~ cVl (13) 

for some c G (0, oo). 

Proof. If the density A of +1 sites is larger than p s c lte , there exists almost surely an 
infinite cluster of +1 sites at time zero. Every doubly-infinite path or closed loop of sites 
of the same sign is stable for the dynamics. By standard percolation arguments, there are 
loops of stable +1 sites that break up the lattice into finite subsets, and the probability 
that such subsets are large is stretched-exponentially small. To be more precise, if D is 
the subset that contains the origin and is surrounded by the smallest such loop, then 

PA\Do\>N)<e- c '^ (14) 

for some d G (0, oo). To sketch the arguments that lead to this, let Rl be any L x 3L 
rectangle in M 2 and consider the event A(Rl) that the set Rl(\M contains a plus-crossing 
joining its two sides of length L. By standard arguments ( 1)15] ; see also |28j and |2J 



P a o{A{R L ))>l-e^ L (15) 
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for L large enough and some (3 G (0, oo). Let B(L) be the set ([-L/2, L/2] x [-L/2, L/2])f) 
HI and S(L) be the event that there is a plus loop in B(3L) \B(L) surrounding the origin. 
Then, writing B(3L) \ B(L) as a union of four L x 3L rectangles and using (|T5|), we have 

PAW > (3L) 2 ) < l-P a o(S(L)) 

< e~ c " L (16) 
for L large enough and some c" G (0, oo), which yields (14). 



By Proposition |2]3|, the maximum number of spin flips allowed inside D is 8 |D |, 
where \D \ is the cardinality of D Q . Moreover, as long as a single unstable site is present, 
it will flip at its next clock ring. It follows that an upper bound for the time for the sites 
in Dq to fixate is given by a sum of 8 \Dq\ independent exponential (mean one) random 
variables T%, .... More precisely, for any a > 0, 

P{t) < Po°(\D \ > [at]/ 8) + max Prob(Ti + ...+T n >t) 

l<n<[a t] 

= P o0 (\D Q \>[at]/8) + Pi6b(T 1 + ... + T [at] >t). (17) 

By choosing < a < 1, the first term at the end of ([17]) can be seen to be exponentially 
small in \rt, using (|14|), and the second term is even smaller (exponentially small in t) by 
standard large deviation arguments. □ 

The next result is a slight improvement of the previous theorem. 

Theorem 2. There exists A < p s c lte such that i/A 6 (A, 1] in the homogeneous ferromag- 
net, then 

p(t) < e~ cVi (18) 

for some c G (0, 00). 

Proof. To prove the theorem, we will use a general result of Aizenman and Grimmett 



P|, [15j on enhanced percolation. To do this, first of all let us partition the hexagonal 
lattice into two sublattices A and B (with the origin of HI in B) in such a way that all 
three neighbors of any site in A (resp., B) are in B (resp., A). Any two sites of either A or 
of B have no edge of the original hexagonal lattice in common. By joining two sites of A 
whenever they are next-nearest neighbors in the hexagonal lattice (two steps away from 
each other), we get a triangular lattice (the same with B). For each site x G A, write the 
exponential random variable representing the time of the first ring of the Poisson clock at 
x as a sum of three independent identically distributed random variables (this is possible 
because the exponential distribution is infinitely divisible) . Assign each one of these three 
random variables to one of the three neighbors of site x (such neighbors belong to the 
sub lattice B). We now construct our enhanced percolation process in the following way: 

• A site in the sublattice A is open if its spin is plus. 
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• A site in the sublattice B is open if 



1. its spin is plus, or 

2. the random variable representing the first clock ring of that site is smaller than 
each of the three random variables it was assigned from its three neighbors in 
A and smaller than one. 

This last condition implies that that site in B will make its first attempt at flipping before 
any of its three neighbors in A. The enhancement only takes place at sites of the sublattice 
B and is activated at site x with strictly positive probability. 

The enhancement is easily seen to be essential, as defined in [Q, and therefore the 
result of Aizenman and Grimmett [|TJ (see also Section 3.3 of ]T3|) can be applied. Thus, 



the critical value for the density of +1 sites in cr° to yield percolation of open sites 
is strictly lower than that for independent site percolation: p e ™ h < pf te . To conclude the 
proof we need to show that A > p e ™ h implies stretched exponential decay of pit)- This is 
so because we claim that any loop of open sites in EI is a loop of plus sites at time one, 
from which stretched-exponential decay with exponent 1/2 follows by using essentially 
the same arguments already used in the proof of Theorem |I[ 

To see why any loop of open sites at time zero becomes a loop of plus sites at time 
one, notice that by time one each open site of the lattice B which was not plus at time 
zero will have attempted a spin flip before any of its neighbors and therefore will have 
flipped to plus. □ 

Remark 2.4. We conjecture that p(t) decays exponentially for all values of the initial 
density of +1 spins X, including 1/2 (see also fcU\ , \23[]). This is also supported by simu- 
lation results of Howard / {T^| , [7^ /. 



3 Synchronous dynamics 

In this section we study a different kind of dynamics, in discrete time, for the homogeneous 
ferromagnet on EL This is the zero-temperature case of a dynamics studied by Domany 
The rules (but not the timing) for updating the spins are the same as in the previous 
section and we will use the same notation. The hexagonal lattice is partitioned into two 
sublattices A and B (as in the proof of Theorem H) in such a way that the set M{x) of 
all three neighbors of a site x in A (resp., B) is in B (resp., A). By joining two sites of A 
whenever they are next-nearest neighbors in the hexagonal lattice (two steps away from 
each other), we get a triangular lattice (the same with B). The synchronous dynamics is 
such that all the sites in the sublattice A (resp., B) are updated simultaneously. Since 
this is a discrete time dynamics, our stochastic process will be denoted by a n , with 
n G {0, 1, 2, . . . }; cr° is still chosen from the Bernoulli distribution P a o with density A of 
+ 1 spins. The first update a 1 will be for the sublattice A. 
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Proposition 3.1. If A > 1/2 (resp., < 1/2), there is percolation of +1 (resp., —1) spins 
in a n for any n G [1, oo] . 

Proof. We will only give the proof for A > 1/2, since the case A < 1/2 is the same 
by symmetry. If A > 1/2, since the critical value for independent (Bernoulli) percolation 
on the triangular lattice is exactly 1/2, there is at time zero percolation of +1 spins in 
the two triangular sublattices A and B. Any site x G A with two +1 neighbors in B 
will become +1 when it updates, but since +1 spins percolate (and form doubly-infinite 
paths) in B, at time 0, when the sites of A are updated at time 1, doubly-infinite paths 
of +1 spins will be created in the hexagonal lattice and these are then stable. Therefore, 
for n > 1, there is percolation of +1 spins in the hexagonal lattice. □ 

Theorem 3. For any \, a n converges to a°° exponentially fast for the synchronous dy- 
namics in the sense that 



for some c G (0, oo), where p A {n) is the probability that a deterministic site in A flips 
after time n and similarly for p B {n) . 

Proof. Without loss of generality, we assume that < A < 1. First of all we make 
the following observations: 

• the values of the spins in the sublattice A at time zero are irrelevant, since at time 
1, after the first update, those values are uniquely determined by the values of the 
spins in the sublattice B, 

• once the initial spin configuration in the sublattice B is chosen, the dynamics is 
completely deterministic. 

We now concentrate on the sublattice B. It is easy to see that the deterministic 
dynamics in B (observed at even times 2n) is of nearest neighbor type. We also claim 
that each site x G B can flip at most one time. To see this, assume without loss of 
generality that a —1 spin at x G B flips to +1 at time 2n with n > 1. For this to happen, 
two of the three H-neighbors of x (yi,V2 G A) must be +1 at time 2n. Therefore, at time 
2n — 1, yi and y<i each need to have two +1 H-neighbors in B. This implies that at time 
2n — 1 all the sites in the hexagon containing x,y±,y2, except x, are +1. When, at time 
2n, the spin at x flips to +1, a stable loop is formed, and site x fixates. 

Let us consider a loop 7 in the triangular sublattice B, written as an ordered sequence 
of sites (2/0, 2/1, ... , y n ) with n > 3, which are distinct except that y n = y . For i — 1, . . . , n, 
let Q be the unique site in A that is an H- neighbor of both and yi. We call 7 an s-loop 
if £1, ... , ( n are all distinct. Similarly, a (site-self avoiding) path (yo, yi, . . . , y n ) in B is an 
s-path if Ci, • • • , (n are all distinct. Notice that any path in B (seen as a collection of sites) 
contains an s-path. We will abuse our terminology slightly and also call a doubly-infinite 
s-path an s-loop. An s-loop of constant sign is stable for the dynamics since at the next 




(19) 
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update of A the presence of the constant sign s-loop in B will produce a stable loop of that 
sign in the hexagonal lattice. A triangular loop x±, x 2 , £3 G B with a common H-neighbor 
C G ^4 is called a star; it is not an s-loop. A triangular loop in B that is not a star is 
an s-loop and will be called an antistar, while any loop in B that contains more than 
three sites contains an s-loop. A connected cluster in B that does not contain any loops 
other than stars (i.e., it doesn't contain s- loops) will be called an s-tree (notice that this 
does not correspond to the usual definition of a tree). An s-tree can contain more than 
one star, but any two distinct stars cannot have sites in common or they would together 
contain an s-loop. 

Consider a specified site x (e.g., the origin) in the triangular sublattice B and denote 
by its £>-cluster of constant sign (at some time). The site x can be of three types: 

• x belongs to an s-loop or to an s-barbell which consists of two s-loops of constant 
sign connected by a path (and therefore by an s-path) of the same sign, in which 
case it is called frozen, 

• x has only one B- neighbor of its sign or else exactly two B- neighbors y\,y 2 of its 
sign and x, yi, y 2 share a common H-neighbor (so that (x, yi, y 2 ) is a star), in which 
case it is called hot, 

• x is neither frozen nor hot, in which case it is called warm. 
Notice that: 

• a frozen site will never flip at any later time; a hot site will flip at the next time 
step; a warm site will not flip at the next attempt, but might flip at a subsequent 
time, 

• any cluster that is not completely frozen must contain at least one hot site (because 
the removal of all frozen sites results in one or more s-trees), 

• a cluster can grow but every newly added site will be frozen (as shown in the 
argument above that each site can only flip once), each newly added site being 
previously a hot site of a cluster of opposite sign. 

To analyze when a site x G B will flip for the last time, we consider its £>-cluster at 
time zero. According to the above considerations, if x is hot at time zero, then it will 
fixate at the first update of £>; if it has four or more £>-neighbors of its sign at time zero, 
then it belongs to an s-loop and is frozen. There are only three possibilities left to be 
considered: 

1. x has exactly two ^-neighbors, y\ and y 2 , of its own sign and such that x,y±,y 2 
don't form a star. 

2. x has exactly three B- neighbors, y±, y 2 and 2/3, of its own sign and such that none 
of the combinations x, y^ yj, for i, j = 1, 2, 3, forms a star. 
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3. x has exactly three £>-neighbors, y\, y-z and y%, of its own sign and such that x, y±, y2 
form a star and y% is not a i3-neighbor of yi or y 2 . 

Now, let U be a specific nonempty subset of the set Af B (x) of all six £>-neighbors of 
x. We define the "partial cluster" C® v to be the set of sites y G B such that there is a 
(site-self avoiding) path xq = x, xi, . . . , x n = y of constant spin value (at time zero) in 
B with x\ G U; i.e. is the set of sites in C® that can be reached from x by paths 

that start off by going from x to a neighbor in U. Note that if we are considering one of 
the cases 1, 2 or 3 described above, we can define a branch of C® as C% v with U either 
a singleton or, in case 3, also the doubleton {2/1,2/2}- In case 1, C® has two branches; in 
case 2, it has three branches; in case 3, it has two branches, one of which has a doubleton 
U. For x not to be frozen, the branches of 0~ must be distinct (otherwise x would be 
in an s-loop), and all but one of them must be s-trees, i.e., must contain no loops other 
than stars (otherwise x would be in an s-barbell). In this case, x may eventually flip, and 
the time at which it will do so is bounded by the length of the longest s-path contained 
in one of the branches of C® that is an s-tree at time zero. 

To complete the proof, it suffices to show that, for any fixed U, there is some (3 > 
and K < 00 such that 

P(\C%v\ ^ n and C xu contains no antistar) < Ke~ Pn . (20) 



To prove (£(]), we partition B into disjoint antistars and denote by r the collection of 



these antistars. We do an algorithmic construction of C® v (as in, e.g., [0), where the 
order of checking the sign of sites is such that when the first site in an antistar from r is 
checked (and found to have the same sign as x), then the other two sites of that antistar 
are checked next. Without loss of generality, we assume that cr° = +1. Then standard 
arguments show that the probability in (|20|) is bounded by K (1 — A 3 /™/ 3 -*, which yields 
(H) for p B {n). 

If a site z belongs to the triangular sublattice A, it is clear that once its three H- 
neighbors have fixated, z will certainly be fixated at the next update of A, and therefore 
the exponential bound (|19|) holds also for p A {n). □ 

Theorem 4. If A = 1/2 in the synchronous dynamics : 

1. For n G [0, 00], there is no percolation in o n of either +1 or —1 spins, for almost 
every a . 

2. There is infinite mean cluster size in o n for any n G [1, 00] : for any x G H, 

£7(10,(71)1) = 00. (21) 



Proof. The proof of the first claim is the same as the proof of Proposition |2.1| . To 
prove the second claim, let us first notice that C x {n) n B = C®{n) for any x G B and odd 
n > 1 and C x (n) H A = C£(n) for any x G A and even n > 2. Thus percolation in B for 
odd n > 1 (resp., in A for even n > 2) would imply percolation in the hexagonal lattice 
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at the same n. Therefore, by the first part of the theorem, at any odd n > 1 (resp., even 
n > 2) there is no percolation in B (resp., A). By a theorem of Russo 



(see also |27| ) 

applied to the triangular lattice, this implies that the mean cluster size of the +1 and — 1 
clusters in, say, B at odd n > 1 diverges. It follows that, for any odd n > 1 and x G B 



E{\C x {n)\) > E{\C B x {n)\) 



oo. 



Now, if |C x (n)| 7^ 1, then C x {n) = C y (n) for some y G Af(x) and so 

\C x (n)\< \Cy(n)\- 

yt=M(x) 



(22) 



(23) 



Taking expectations and noting that each neighbor y of x G B is in A, we see that for 
odd n > 1, ( 22|) for x G B implies (^l|) for all x G H. The proof of ([H]) for even n > 2 is 
similar with ^4 and i3 interchanging their roles. □ 



4 Disordered ferromagnet 

In this section we study a disordered ferromagnet with i.i.d. couplings J x<y uniformly 
distributed between and 1. We will also have an external magnetic field h > 0, so that 
the (formal) Hamiltonian is 

U = - 7^ Jx^CTxCTy ~ h Ox- (24) 

The initial density of +1 sites is A G [0, 1]. Let P = P a o x x Pj denote the joint 
distribution on the space of cr°'s, tu's and J7"'s. 

Remark 4.1. T7ie theorems and proofs of this section are valid for more general choices 
of the common distribution of the J XtV 's than uniform on [0, 1]. For example, Proposition 
is valid for any distribution on [0, 1) and Theorem [| is valid for any continuous 
distribution on (0, 1) that is symmetric about 1/2. 

Proposition 4.1. In the context of the disordered ferromagnet just described, if h > 2 
and X G [0,1], almost surely, lim^oo a x = +1 for all x. Moreover, fixation happens 
exponentially fast in the following sense: 

p(t) < e- ct (25) 

for some c G (0, oo). 
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Proof. Let p(t) denote the P a o xP w x Pj probability that the origin is —1 at any time 
after t. We claim that to prove both conclusions of the proposition, it suffices to show 
that p(t) < e~ ct . To see this, first note that p(t) — > implies that the origin is eventually 
+ 1 almost surely, and so by translation invariance, the same is true for any site; then 
note that p(t) < p{t) (since if a site flips then it must be minus either just before or just 
after the flip). 

Now notice that there is a positive density p of sites x such that J2 y eAf(x) J%,y < ^- 
If the origin is one of those sites, then either it is +1 and will remain such, or it is — 1 
and will flip to +1 the first time its clock rings, which happens exponentially fast. If the 
origin is not one of those sites, there is a closest site X such that J2 y eAf(x) Jx,v < h an d> 
since p > 0, 

Pj(\\X\\>n)<e- c ' n (26) 

for some d G (0, oo), where denotes the number of steps between the origin and 

site X along a shortest path. Next, notice that if a site z has at least one neighbor y 
with o y = +1 when the clock at z rings, a z will flip to +1 if it is —1 and will remain +1 
otherwise. Then we have, for any a > 0, 

p(t) < Pj(\\X\\ > [at]) + max Prob(T + T x + . . . + T n > t) 

0<n<[at] 

< P J (\\X\\>[at]) + Prob(T + T 1 + ...+T [at] >t) (27) 

where T ,T 1; . . . are independent exponential (mean one) random variables. By choosing 
a small, both terms above can be seen to be exponentially small in t, the first one using 
(p6|), the second one by standard large deviation arguments. □ 

When < h < 2, it is easy to see that a°° = lim^oo a 1 exists almost surely, but is not 
constant, that is, for any A G (0, 1) there is a positive density in cr°° of both +1 and —1 
sites. In this case, we can prove (in the next theorem) exponential fixation (in the above 
sense) only for h > 1.5, although we conjecture that exponential fixation happens for all 
values of h down to h = 0. 

Theorem 5. In the disordered ferromagnet with A G [0, 1], if h > 1.5, 

pit) < e- ct (28) 

for some c G (0, oo). 

To prove Theorem [|, we need the following definitions and lemmas: 

Definition 4.1. For each site x, let t(x) = {r n (x) : n G N} denote the arrival times 
(i.e., times of clock rings) of the Poisson clock associated with site x, arranged so that 
Tk(x) < Tk+i(x) for each k. Let r + (x,t) = inf{rk(x) : Tk(x) > t] be the next clock ring at 
x after time t. 
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Definition 4.2. WE We call a sequence of sites (x\, . . . ,x n ) with Xj G M{xj-\) for 
j = 2, ... ,n a plus-cascade (of length n) if, for a sequence of times ti, . . . ,t n defined by 
t\ = Ti(xi) and, for k > 1, tj. = T + (xk,tk-i), we have: a Xk flips from — 1 to +1 at time 
t = tk and, for k < n, does not flip again until possibly after time tk+i, and for k > 1, 

tk > n(x k ). 

Similarly, we define minus-cascades. A sequence of neighboring sites (x\, . . . , x n ) is a 
cascade if it is either a plus- or a minus-cascade. If site x does not flip at ti(x) (x's first 
clock ring), then no cascade begins at x and we say that x has an empty cascade. Notice 
that in the disordered models, cascades can split and merge. As an immediate consequence 
of the definition, a plus-cascade (x\, . . . ,x n ) has the two following properties, which are 
useful to keep in mind: 

• for 1 < k < n, Xk and $k-\-l £11*6 nearest neighbors, 

• for 1 < k < n, at the time t — tk when a x . flips to +1, crl = +1 and o~i = —1. 
Similar properties are valid for a minus-cascade. 

Lemma 4.1. Every flip belongs to some cascade. 



Proof. Consider a spin flip from (say) —1 to +1 at site x at time t. If t = Ti(x), 
then the flip we are considering belongs to a cascade starting at x. If t > Tx(x), then 

t = Tk(x) for some k > 2. In that case <r* + = — al (k ~ 1) x + , which implies that at least 
one site y± G N(x) experienced a spin flip from —1 to +1 at some time t\ G (r(fc_i)(x),t) 
and then did not flip again before time t. We repeat the same procedure for y% and 
construct inductively a sequence of sites (yo = x, y\, . . . , y n ) and a sequence of times 
to — t > t\ > ■ ■ ■ > t n . The procedure stops when t n = Ti(y n ). It should be clear 
from the construction that (y n , . . . , y x , y = x) is a cascade for the sequence of times 
t n , . . . , t\, to = t. □ 



Lemma 4.2. In the setting of Theorem each site can flip at most 7 times. 

Proof. By the translation invariance of the model, it is enough to prove the claim for 
the origin 0. In order to do that, we will show that a minus-cascade can have at most 
length 2. Consider a minus-cascade whose first two sites are respectively xo and yo- Call 
x\ and X2 the other two neighboring sites of xq, and y\ and yi those of yo- When xo 
flips for the first time, from +1 to —1, by the definition of a minus-cascade, crll^ = +1, 
because yo is the second site of the cascade. Then, because of the strength of the field 
(h > 1.5), for a Xo to flip from +1 to —1 at time t\ = Ti(xo), the two following conditions 
are necessary: al 1 } 10 ^ = al\ — — 1, and J X(hyo < 0.5. Now, for a yo to flip from +1 to 
— 1 at time t 2 ? since J xo ,y < 0.5, it must be the case that o~ yi = o~ y2 = — 1. Therefore, 
the cascade cannot proceed further, since when a yo flips at time t 2 ? all its neighbors are 
already —1. 
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Consider now the origin 0. It can flip from +1 to —1 at most three times, due to a cas- 
cade starting at the origin itself plus at most two cascades starting at two of its neighbors 
or else due to cascades starting at its three neighboring sites (if there is a minus-cascade 
starting at the origin and then o"o flips back to +1, that flip belongs to a plus-cascade 
involving some y G A/"(0) and y cannot then be the start of a minus-cascade reaching 0). 
In fact, a minus-cascade starting further away would not reach the origin. Therefore, the 
origin can flip at most 7 times. □ 

Proof of Theorem |5|. Once again, we partition the hexagonal lattice into the two 
triangular sublattices A and B as in the proof of Theorem |2|. In addition to the Poisson 
clock C x of rate 1 at site x G B, we also assign three more clocks C XiV with rate 1/3 
each and associated respectively with the neighboring sites y G A/"(x), all of them being 
in A. (We remark that these clocks C XjJ/ are different from the ones used in the proof of 
Theorem ^|.) All the clocks C x , C Xiy for x G B and y G Af(x) are independent. We now 
define a new dynamics according to the following rules: 

• A site x G B attempts a spin flip when C x rings. 

• A site y G A attempts a spin flip whenever any one of the clocks C XtV assigned to 
one of its neighbors x G J\f(y) (in sublattice B) rings. 

• The rules for accepting a spin flip are the same as before. 

Since the clocks C x ,y used for sites y in A have each rate 1/3, the resulting dynamics 
is the same as before. We fix a time T > and say that x G B is a good site if the two 
following conditions hold: 

L T, y eX(x) J *,y < K and 

2. during the time interval [0, T], C x rings and following that all three clocks C x>y for 
y G Af(x) also ring during [0,T]. 

A site that is not good will be called bad. Clearly, because of condition 1 and the external 
magnetic field, good sites that are ever +1 are then stable for the dynamics. Condition 
2 implies that a good site x G B that is —1 at time zero will have flipped to become a 
stable +1 by time T, and moreover it implies that a site y G A with two good sites as 
H-neighbors will be a stable +1 by time T, regardless of its value at time zero. 

The event corresponding to condition 1 has probability strictly larger than 1/2 be- 
cause of the distribution of the couplings and our assumption that h > 1.5. The event 
corresponding to condition 2 has probability approaching one as T — > oo. So, choosing T 
large enough, we will have that the probability that a site x is good is strictly larger than 
1/2. Thus good sites percolate in the triangular lattice B and therefore, for x G B (say, 
the origin), if B® is a £>-cluster of bad sites, 

(P u x Pj){\bH\ > n) < e-* n (29) 
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for some c" G (0, oo). Moreover, the boundary dB® of the cluster B® is the union of (one or 
more) £>-connected loops of good sites. Let d£ = {y G A : 3 Ci; C2 £ 9B® s.t. | \d — y\ \ = 
IIC2 — y\\ — 1} an d define d x = dB® U 9^ . Then d x is the union of (one or more) El- 
connected sets and it completely surrounds the H-connected set B x = B® U B£, where 
B* = {y G -4 : 3Ci,C 2 e Bj? s.t.||Ci -y|| = ||Ca — 3/|| = 1}- Clearly |£y < 3 |Bf | and 
therefore 

(P w x Pj)(\B x \ >n)< e~ c "™/ 3 = e" c '" (30) 

for some d G (0, 00). 

Because of condition 2 above, by time T all the sites in d x have fixated: the sites 
in dB® because they are good sites, those in d£ because they have two good sites as 
H-neighbors. To conclude the proof, observe that if the origin has not fixated by time T, 
it must belong to a connected set B surrounded by fixated sites and with 

(P UJ xP J )(\B \>n)<e- c ' n . (31) 

By Lemma [4.2| , the maximum number of spin flips allowed inside Bq is 7 \Bq\. Moreover, 
as long as a single unstable site is present, it will flip at its next clock ring. It follows 
that an upper bound for the time for the sites in B to fixate is given by a sum of 7 |-So| 
exponential (mean one) random variables Ti,T 2 , .... More precisely, for any a > and 
t > T, 

P(t) < (P u xPj)(\B \>[a(t-T)}/7) 

+ max Prob(Ti + . . . + T n > (t-T)) 

l<n<[a(t-T)] 

= (P UJ xP J )(\B \>[a(t-T)]/7) 

+Prob(T! + . . . + T [a(t _ T)] > (t - T)). (32) 

By choosing < a < 1, both terms in (|32|) can be seen to be exponentially small in t, the 
first one using (pi|) , the second one by standard large deviation arguments. □ 
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